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Abstract 


We show an alternative proof for the existence of weak solutions to equations describing turbulent flows of fluids. 
The proof proposed by one of the authors in a previous paper (cf. [A. Swierczewska, Large Eddy Simulation. 
Existence of Stationary Solutions to a Dynamical Model (submitted for publication). Preprint TU-Darmstadt 
no. 2314, http://wwwbib.mathematik.tu-darmstadt.de/Math-Net/Preprints/Listen/shadow/pp2314.html]) based on 
more classical methods. We will use Young measures, which allow us to shorten significantly the limiting procedure 
in the nonlinear terms and generalize the statement. 
© 2005 Elsevier Ltd. All rights reserved. 
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1. Introduction 


In the paper [1] the author proved the existence of stationary weak solutions to the equations 
describing turbulent flows in the three-dimensional torus T?: 


v: Vu — div A (y, Vv) — vAv+ Vq =f, 


1.1 
div v = 0, (1.1) 
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where v : T? — R? is the velocity, q : T? — R is the pressure and V°v = (Vu + Vv") denotes the 
symmetric part of the gradient. The operator A is a nonlocal operator, continuous w.r.t. to both variables 
y = (0, VŽ Ù, vv, |V*u|V5v) and V*v. The tilde denotes a convolution with some smooth function ø, 
namely a(x) = fgs uY) — y) dy. 

The equations come from the Large Eddy Simulation (LES) technique for incompressible flows. The 
idea of LES has its origin in numerics. Typical for turbulent flows are very different scales, which lead to 
an increase in the number of numerical operations needed to compute the solution. The LES technique is 
based on choosing the scales for which the exact solution is computed directly and the scales for which 
the solution is modelled. The selection of scales is obtained through filtering, i.e. convolving with some 
function (filter) the Navier-Stokes equations. Adding a constitutive relation representing the contribution 
of small scales into the flow leads to the Eq. (1.1). The modelling is clearly described in [2—4]. The 
monograph [4] is a wide description of different modelling techniques in LES. For mathematical results 
for one of the other approaches we refer to [5,6]. 

We recall the existence theorem and briefly present the idea of the proof. The set of smooth functions 
on the torus T? can be identified with the set of periodic smooth functions with some period L € (0, 00). 
Therefore 2 = (0, L)? is a cube of side L in R3. By cR’) we denote the set of C% (R3) functions, 
which are periodic with period L > 0 in each ith direction, i.e., u(x + Le;) = u(x), i = 1, 2,3. In the 
following V denotes the closure of the space 


v= fu : u € CH, R’), divu =0, f wax =o} 


with respect to the norm ||u||y = (So |Vu]? dx)3. Its dual space will be denoted by V’. 


Theorem 1.1 (Existence). Given f € V there exists a weak solution to the stationary problem (1.1), 
i.e. the equation 


[9-00 + AG, VD: 19+ Vv: Vo- fodar =0 (1.2) 
R 


is satisfied for all p € V, where y = (ù, VSÙ, vv, [V5 u]V5v). 


Most of the difficulties in the existence proof are concentrated in passing to the limit in the nonlinear 
term A, given by 


A(y, V°v) = c)|V v| Vā v, 

with the function c continuous with respect to all variables and satisfying the condition 
O<a<c<B<o. 

With the help of the properties of convolutions, the strong convergence of the subsequence 
y™ —> yin L*(Q) 

is obtained, where y” = (v0", V°u", v@v", |V5u"|V5v"), y = (Č, VS, Gv, Z) and |V v" |V v" —> x 


in L3(Q). First, with the help of the Minty Browder trick, a better characterization of the weak limit 
c(y"*) [Vou [Vou — e(y)| V v| V5 v in L3(2) is achieved. 
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Noting additionally that the structure of Galerkin approximations yields that 
im f eV ax = f eGyiv'ur ax. 
n—> Oo Q Q 


allows us to show the convergence of appropriate weighted norms and hence the strong convergence of 
gradients, which completes the proof. 

In the following, we show an alternative way of passing to this limit. It is much shorter; however, 
it uses more advanced techniques. Below we formulate a more abstract theorem, which solves also the 
problem of passing to the limit in the turbulent term of system (1.1). Moreover this approach allows one 
to solve the problem in bounded domains, which is more complicated because of convolving near the 
boundary. For the sake of conciseness we do not present this case in detail. 


Theorem 1.2. Let OQ C R” be a measurable set of finite measure and let an operator A(x,s,&) : 
NR x R” x R” — R” satisfy the following conditions: 


(i) A(x, s, €) is a Carathéodory function (measurable w.r.t. x and continuous w.r.t. (s, E)). 
(ii) Forallx € 2,5 € R” and £1, & € R",& 4 & 
[A@, s, &1) — A, s, &2)] - [&1 — &2] > 0. 
(iii) There exist positive constants c1, c2 such that for p > 1 it holds that 
A(x, s, E) -E > cill” 
and 


|A(x, s, Ẹ)| < c2l£|P7t. 


Let y” : 22 > R” and z” : Q — R” be sequences of measurable functions such that 


(iv) y" > y a.e. in Q, 
(v) z” — z in LP (Q) and A(x, y”, z”) = Ain LPT (0), 
(vi) 


tim sup f A(x, y”, z”): z” dx <f A- zdx. 
R R 


n—> 00 


Then there exists a subsequence of (z") such that 
z” > zin L’ (Q). 


Remark. In the proof of Theorem 1.1 in [1], the Minty Browder trick does not give a 
good characterization of the term x because the operator A(y, V*°v) is not pseudomonotone. 
However Theorem 1.2 shows that A is of type M (i.e. if X is a reflexive Banach space, A : X > 
X’, u” — u, Au” — x, limsup,_,,,.(Au",u") < (x,u), then Au = x). For more details on the theory 
of pseudomonotone operators we refer the reader to [7]. The idea of Theorem 1.2 is influenced by the 
compactness properties of Leray—Lions operator (cf. [8, Lemma 5, p. 190]). However their proof follows 
quite a different path to ours. 
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2. Young measures tools 


For the convenience of the reader we collect below all the necessary tools concerning Young measures 
used in the proof of Theorem 1.2. For more details and the proofs, we refer the reader to [9, Corollaries 
3.2-3.4]; see also [10,11]. 


Lemma 2.1. Suppose that the sequence of maps zi : Q > R? generates the Young measure v. Let 
F : R x R? —> R be a Carathéodory function (i.e. measurable in the first argument and continuous 
in the second). Let us also assume that the negative part F~ (x, zÍ (x)) is weakly relatively compact in 
L! (Q). Then 


liminf | Fados ff F(x, A) dv, (à) dx. 
R N JRI 


J> œ 


If, in addition, the sequence of functions x +> |F|\(x, z/ (x)) is weakly relatively compact in L! (2), then 
FO, z/()) > J F(x, A) dv (A) in L'(Q) 
Rd 


Remark. The second part of the above theorem can be easily extended to vector valued functions F. 


Lemma 2.2. Let u/ : N > RI, vi : Q > RY be measurable and suppose that u —> u a.e. while v/ 
. . [A 

generates the Young measure v. Then the sequence of pairs (u! , v?) : R —> RI+4 generates the Young 

measure x +> Sy(x) 8 Vx. 


Lemma 2.3. Suppose that a sequence z/ of measurable functions from N to R? generates the Young 
measure v : N > M(R¢). Then 


z! — z in measure if and only if vy = ôz(x) a.e. 


3. Proof of Theorem 1.2 


We apply Lemma 2.1 to the function A(x, y”, z”) - z”. The coercivity condition from assumption (iii) 
of the theorem assures that the negative part of this function is equal to zero; thus it is certainly weakly 
relatively compact in L! (2). This allows us to conclude that 


imsup f ACs, »",2")-2" ds =f f M8 8) E dials, E) dx (3.3) 


n—> 0 


where uy is the Young measure generated by the sequence (y”, z”). However according to Lemma 2.2, 
we are able to characterize this Young measure more precisely. Since the sequence y” — y a.e. and z” 
generates the Young measure vy, the Young measure uy generated by this pair satisfies wy = ôṣ5(x) @ Vx- 
Therefore we can integrate 


J] A@s E) Eda Dda ff EEEE T EO (3.4) 
R Jem xR 2 JR" 


In the same way we obtain 


J] aœ s Dans Da = ff TDA de. (3.5) 
2 JR™ xR" RJR" 
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Since the sequence |A(x, y”, z”)| is bounded in L?-T(Q), it is weakly relatively compact in L'(), 
which implies A = Ji A(x, s, €)dux(s, €). Thus, from (3.3)-(3.5) and assumption (vi), the following 
inequality holds: 


[ Í. ATE: I. £ dv, (Ẹ) dx > Lf. A(x, 5(x), E) £ dv (&) dx. (3.6) 
From the monotonicity of A w.r.t. the last variable we can deduce that 
j) / h(x, &) dv,(&) dx = 0, (3.7) 
Q R” 


where h is defined by 


h(x, §) = [Aes 50.8) -A (3%. fmo) i E = [sone] 


Simple calculations imply that 


Jf nox, sydv(gax = ff Ax, F(x), &) « & dv, (£) dx 
N JR” N JR" 


- | | Acse.aancy- f Eed, 

‘2: R” R” 

which, together with (3.6), assures that 
Jf h(x, £) dvx (Ẹ) dx < 0. (3.8) 
Q JR" 


Then, (3.7) and (3.8) imply that = h(x, €) dv, (£) = 0 for a.a. x € 2. Moreover, since vy > 0, we have 


wpet fS gao). 


Note that the single point in the right-hand side set is located a.e. in the point z(x), where z is the weak 
limit of the sequence (z”). Finally we can conclude that vy = ôz(x) a.e. A direct application of Lemma 2.3 
yields that z” — z in measure. Then there exists a subsequence of (z”) such that z” — z a.e. Using the 
information that v, = ôz(x) together with Lemma 2.1 and assumption (vi) yields 


tim sup f A(x, y",z")-z" dx <f A(x, y,z)-zdx < timinf f A(x, y”, z”) z” dx. 
Q Q n—> oo Q 


n—> 00 


Hence we can set g” = A(x, y”, z”) - z”, g = A(x, y, z) - z and claim that 


g” >00, g ELA), | a> f sar, g” > ga.e.in R. 
R R 


Noticing that 


J ie" — gldr = fe" -g)dx +2 f max{(g — g8”), 0} dx 
R 2 a 


we conclude by Lebesgue’s Dominated Convergence Theorem that A(x, y”, z”)z” — A(x, y,z)z in 
L'(Q). Thus, by Vitali’s Theorem, it is uniformly integrable. Due to the coercivity condition, also the 
sequence |z”|? is uniformly integrable. Using again Vitali’s Theorem yields that z” — z in L?” (Q), 
which completes the proof. 
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4. New proof of Theorem 1.1 


We recall the theorem concerning the existence of a solution to the Galerkin approximation: 


Theorem 4.1. For any given f € V' there exists a solution à}, ..., àn (and therefore v” € V" = 
span{@1,..., @n}) to the approximate problem 
n 
v" (x) = So Mar, ER, (4.9) 
r=1 
b(v", v”, wr) +f A(y", Vv") - Va, dx + v(Vv", Vo) = (f, or) (4.10) 
Q 
dui ,,,. 

forr =1,...,n, where we denote b(u, v, w) = fo Wj ae Wi dx. 


From the energy estimate we conclude the convergence (even if the convergence holds only for a 
subsequence, throughout this proof we will not use subsequence indexes) 


v” > vin V, 
which implies 
v” > vin L4(2) for 1 <q < œ, 
b(v”, v”, ġ) = b(v, v, œ) forall@ € V 
and we can conclude the existence of A € L? (2), y € LO(2) such that 
Aly", Vu") > Ain L2(Q), y” > Fin L(A) 


and we recall the notation y” = (v", Vy", unun, IVs v” |V v”), y = (0, V*0, vv, x). To achieve the last 
convergence we use the regularizing properties of convolutions. Hence we apply now Theorem 1.2 to 
the operator 


A(y, Vv) : (R? x S? x S? x S$) x S? > S3, 


where by S? we denote the set of symmetric matrices in R33. We identify S? with R and (R? x S? x 
S? x S?) x S? with R?! x R6. In the notation of Theorem 1.2 y” is as denoted above and converges to 
y a.e.; z” = V*u". The operator A does not depend on x and is continuous w.r.t. all other variables; the 
assumptions (1i)—(iii) are certainly satisfied. The assumption (vi) holds while the approximate equation 
yields 


J Arz dx = -viela + f fot ae, 
R R 


and thus due to the lower semicontinuity of the norm it holds that 


tim sup f A(y",z")-z" dx < -vizi + f foax = f A-zdx. 
Q R R 


n—> 0 


Thus according to Theorem 1.2 we obtain that z” — z a.e. in 2 and strongly in L? (9). This improves the 


information on the convergence of y”; namely we can claim that |V°v"|V*u" —> |V*%v|V5vu in LO (2) 


and thus y = y. Finally this yields that A(y”, z”) —> A(y, z) strongly in L3(Q), which concludes the 
proof. 
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